Emergent charge ordering in near half doped Na .46CoO2 
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We have utilized neutron powder diffraction to probe the crystal structure of layered Na I Co02 
near the half doping composition of x =0.46 over the temperature range of 2 to 600K. Our mea- 
surements show evidence of a dynamic transition in the motion of Na-ions at 300K which coincides 
with the onset of a near zero thermal expansion in the in-plane lattice constants. The effect of the 
Na-ordering on the C0O2 layer is reflected in the octahedral distortion of the two crystallograph- 
ically inequivalent Co-sites and is evident even at high temperatures. We find evidence of a weak 
charge separation into stripes of Co +3 ' 5+e and Co +3,5-E , e ~0.06e below Tco=150K. We argue 
that changes in the Na(l)-0 bond lengths observed at the magnetic transition at T m i=88K reflect 
changes in the electronic state of the C0O2 layer. 



I. INTRODUCTION 



The alkali cobaltates NazCo02 have been the subject 
of intense interest as they are a rare example of compet- 
ing interactions on a triangular lattice that can be easily 
tuned by chemical means. Varying the amount of Na (x) 
produces a rich phase diagram which exhibits spin de- 
pendent thermopower (x =0.75) [lj], metal- insulator tran- 
sitions (x =0.5)0, 0], antiferromagnetism and 5K su- 
perconductivity at x =0.3 for a hydrated compound[4] 
. More recently it has been realized both experimentally 
and theoretically that the role of the Na ions goes beyond 
providing a simple means to electronically dope the C0O2 
layer [f| [y, 0]. Rather, the ordering of Na-ions leads to a 
potential that perturbs the C0O2 layer to produce strong 
electronic correlations Q. The role of these correlations 
is still under investigation but it demonstrates that these 
materials can exhibit frustration in two different ways, 
one by the triangular topology of the C0O2 layer and the 
other by the Na induced potential. 

This double frustration is best exhibited at half- 
doping. Here the Na ordering results in a relatively 
simple orthorhombic distortion of the parent hexagonal 
phase in sharp contrast to the complex incommensu- 
rate structures found for higher x compounds [5|. For 
x =0.5 Na-ions order as to form stripes as shown in fig[T] 
while the magnetic susceptibility shows two abrupt de- 
creases (see for example inset in figH]) at T m i=88K and at 
2m2=52K[2, 3] . The first transition is associated with the 
onset of a long range antiferromagnetic ordering^, 0, Q 
while the second transition coincides with a sharp rise 
in the resistivity 0, This second transition has been 
ascribed to be driven by charge ordering (CO) of a t 2g 
electron to form distinct LS Co 3+ (t® g ,S = 0) and LS 
Co 4+ (t2 g ,S — 1/2) ions H. Recent /xSR and neutron 
diffraction measurements [3 0, [T3| propose a magnetic 
structure consistent with this picture, as the magnetic 
lattice comprises of stripes of magnetically inactive Co 3+ 
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and antiferromagnetic (AF) coupled Co 4+ (see fig. [Ha)). 

What is striking in this cobaltate is that the se- 
quence of charge ordering and Neel transitions is re- 
versed (T m 2<T m i), compared to classic charge ordered 
systems such as the manganites[i"l| or magnetite [11] and 
has brought the charge ordering picture into some doubt. 
To reconcile the fact that T„ l2 <T m i, Bobroff et aL [l3j 
propose on the basis of NMR measurements a scenario 
of a successive nesting of the Fermi surface (FS) that 
is coupled to a spin density wave in a way in which 
charge carriers are localized successively with decreas- 
ing temperature. However, this idea has been more re- 
cently disputed as the role of the Na-ions and the crystal 
potential that they impose on the C0O2 layer has been 
theoretically treated better. For example stripe Na or- 
dering induces a weak incipient charge ordering on the 
Co0 2 layer H 0, 0, 

however the mechanism of a pro- 
gressive nesting of the FS with decreasing temperature 
arises directly from the Na-ordering as shown by a Hub- 
bard model using the spatially unrestricted Gutzwillcr 
approximation Q. Here the charge ordering is viewed 
to be driven by the ordering of Na above the antifer- 
romagnetic ordering at T m i=88K [|| llij ]. Nevertheless 
this model successfully predicts the correct antiferromag- 
netic ordering 14j and suggests a charge separation into 
stripes of Co+ 3 ' 5+e and Co +3 - 5 ~ £ , with e ~0.06e@. This 
value represents a very weak charge ordering and is be- 
low the lower limit of detection for charge separation by 
NMRpjl, but in agreement with powder diffraction mea- 
surements that suggest e ~0.12e at 10K[l5l|. 

In this paper we use neutron powder diffraction (NPD) 
over a wide temperature range (2-600K) to probe the 
crystal structure of layered Na:rCo0 2 near the half dop- 
ing composition of x =0.46 over the temperature range 
of 2 to 600K. Our NPD measurements show evidence 
of a dynamic transition in the motion of Na-ions at 
300K which coincides with the onset of a near zero ther- 
mal expansion in the in-plane lattice constants of our 
Nao.46Co02 sample. The effect of the Na-ordering on 
the C0O2 layer is reflected in the octahedral distortion 
of the two crystallographically inequivalent Co-sites and 
is evident even at high temperatures. We find evidence 
of a weak charge separation into stripes of Co +3 5+e and 




FIG. 1: (Color online) Crystal structure of the Nao.4eCo02. 
(a) a projection of the ab— plane showing the crystallograph- 
ically inequivelant Co(l) (dark blue) and Co(2) (light blue) 
sites and the Na(l) (yellow) and Na(2) (orange) sites. The 
Na-atom reside alternatively above and below the C0O2 sheet. 
The Na(l) atom is directly above or below the Co(l) atom. 
The red arrows indicate the idealized magnetic ordering of 
nominal low spin Co 4+ sites for the x =0.5 composition pro- 
posed on the basis of neutron and /^SR measurements |8j. (b) 
A portion of the C0O2 sheet. Note that the Na(2) site resides 
above a triangle defined by the edges of C0O6 octahedra. Se- 
lected Co-0 bond lengths determined at 2K are shown. 



Co +35 ~ £ , e ^0.06e below Tco = 150K, thus confirming a 
more physical sequence of charge ordering and magnetic 
transitions for this compound. We argue that changes 
in the Na(l)-0 bond lengths observed at the magnetic 
transition at T ml =88K reflect changes in the electronic 
state of the C0O2 layer. 

The paper is structured in the following way. In section 
III we discuss the evidence for a weak charge ordering as 
determined from the temperature dependent NPD data. 
The dynamic behavior of Na-ions at high temperature 
and changes in Na-0 bond lengths close to the magnetic 
transition T m \ are discussed in section IV, while the un- 
usual temperature dependence of the lattice constants in 
section V. Discussion and summary are found in section 
VI and VII respectively. 
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FIG. 2: Rietveld refinement of NPD data measured at 1.8K 
from our i=0.46 sample. Here crosses represent the mea- 
sured NPD data, while the continuous line through the points 
represents the calculated diffraction pattern. The difference 
between the data and the model is shown at the bottom of 
the figure. Vertical bars represent expected Bragg reflections 
for the Pnmm structure of this compound. The weighted R- 
factor here is wR p —4.22% and the Bragg R factor 4.3%. In the 
inset we show the magnetic susceptibility of our Nao.4eCo02 
powder sample measured on a SQUID magnetometer under a 
field of IkOe. Arrows indicate two anomalies that have been 
ascribed to T m i=88K and T m2 =52K. 



II. EXPERIMENTAL 

Polycrystalline samples were prepared using standard 
solid state synthesis techniques. The starting stoichiom- 
etry for these samples was Nao.7sCo02. In order to 
deintercalate Na from the lattice to achieve a x ~l/2 
composition, a 5g portion of the x =0.75 sample was 
immersed in a bromine- acetonitrile solution with a 1:1 
Na to Br2 ratio, stirred in solution for 7-14 days and 
washed. The Na/Co ratio of the product was measured 
using neutron activation analysis (NAA) giving a compo- 
sition x =0.46(1). Magnetic susceptibility as a function 
of temperature (xCO) was measured using a Quantum 
design MPMS and was found to be identical to thepub- 
lished literature as shown in the inset of fig. 2 0, 0j- 
Rapid measurements of high resolution neutron powder 
diffraction data were collected from the x=0.46(l) sample 
using the HRPD diffractometer (Ad/d ~ 5 x 10~ 4 ) at the 
ISIS-facility, Rutherford Appleton Laboaratory. Higher 
statistics data suitable for Rietveld refinement were mea- 
sured between 2 to 600K using the high resolution pow- 
der diffractometer E9 (Ad/d - 2 x 10~ 3 , A=1.7973A), 
located at the Berlin Neutron Scattering Center, at the 
Hahn-Meitner-Institut (HMI). Supplementary tempera- 
ture dependent data were also measured from the x =0.75 
sample between 5-300K. All NPD data were analyzed 
using the Rietveld method which allowed us to measure 
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lattice parameters, atomic positions and atomic displace- 
ment parameters as a function of temperature. A typical 
Rietveld refinement of the NPD data is shown in fig. [2l 



III. EVIDENCE OF WEAK CHARGE 
ORDERING 

The validity of the reported orthorhombic Vnmm 
structure for x = \ 0, EH] was tested by considering space 
groups that arise from distortions of the parent hexago- 
nal structure P63 jmmc which are consistent with the 
reported orthorhombic unit cell. The program ISODIS- 
PLACE was used for this purpose 17]. This approach for 
the symmetry analysis resulted in space groups that were 
either centered or primitive mono clinic, in both cases in- 
compatible with the diffraction data. Alternately, us- 
ing as a starting point the reported space group Pnmm, 
and testing for related space groups that are compat- 
ible with the diffraction data resulted in primitive non- 
centrosymmetric solutions such as Pmm2 or Pmm2i. Ri- 
etvelds analysis on the basis of these space groups re- 
sulted in somewhat poorer fits than the reported struc- 
ture. Our best modeling of the NPD data between 2- 
450K were obtained using the Vnmm model 0, EH , pro- 
ducing refinements with wR p of approximately 4% (be- 
tween 2- 450K). Contrary to the observation for higher x 
compositions, we find no evidence of incommensurability 
in the NPD data@. 

The crystallography of the Pnmm structure has been 
published elsewhere [3. [HI], however for clarity we illus- 
trate it in fig. [1] and remind the reader that here there are 
two symmetry in-equivalent sites for Na (labeled Na(l) 
and Na(2)) and two for Co (labeled Co(l) and Co(2)). 
The local symmetry for the Co(l) and Co(2) sites dif- 
fers in that the Na(l)-atom lies directly above or be- 
low the Co(l)-atom while the Na(2)-atom resides above 
a space formed between CoOg octahedra as illustrated in 
ng. [lib). 

The effect of the Na-ion potential correlates with the 
distortion of the CoOg octahedra. This distortion is 
shown in fig. [ljb) where we illustrated the difference 
on the CoCvNa coordination and selected Co-0 bond 
lengths. Here the Na-ion potential results in a distor- 
tion of the Co(l)C>6 octahedron, so at 1.8K the six equal 
Co-0 bonds found in metallic x = 0.75 distort to form 
three long bonds (1.90-1.93 A) and three shorter bonds 
(1.84-1.88 A). In sharp contrast the Co(2)06 octahedron 
is more regular with bond lengths varying between 1.88 
to 1.91 A. This indicates that the Co(l) ions experience 
essentially a different crystalline potential than the Co(2) 
ions, as a direct consequence of the Na ordering. In 
fig. [23a, b) we show the temperature dependence of the 
Co-0 bond lengths for the two Co sites. From these data 
it is evident that the larger distortion of the Co(l)06 
octahedron is maintained from low until high tempera- 
ture, while the spread of bond lengths for the Co(2)06 is 
smaller and relatively temperature invariant. These oc- 
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FIG. 3: (Color online) (a,b) Temperature dependence of Co-O 
bond lengths computed from Rietveld analysis of NPD data 
measured on the E9 diffractometer.(c) Bond valence sums for 
the Co(l) and Co(2) sites (filled circles) as a function of tem- 
perature computed from Co-O bond lengths determined from 
Rietveld analysis of the NPD data. The dashed line through 
the data is the average BVS for the two Co sites. BVS values 
obtained from ref. [l5| are also shown. On the same figure we 
plot the octahedral distortion parameter A for the CoOe oc- 
tahedra centered in the Co(l) and Co(2) sites. Dashed lines 
are guides to the eye. 



tahcdral distortions can be quantified by the parameter 
A, where A = y / gfegEg° Eo5g ! ( Co _ O) is the av- 
erage Co-0 bond length and the summation is done over 
the 6 Co-0 bonds. Here A would be zero for a regular 
octahedron with 6 equivalent Co-0 bonds. We find that 
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the Co(l)06 octahedron is more distorted by a factor of 
3 at high temperatures (see fig. EI C )), compared to the 
octahedron centered on the Co(2) site. With decreasing 
temperature however this difference increases to a factor 
of 5 while the distortion saturates below ~150K. As the 
data indicate on the same figure, the distortion of the 
Co(2) octahedron is much less sensitive to temperature. 

At first sight, comparison of these bond lengths at 1.8K 
to the ideal LS Co 3+ and Co 4+ -0 bond lengths of 1.93 
and 1.83 A respectively would suggest that there is no ev- 
idence for integer charge separation between Co(l) and 
Co(2) sites. The analysis of the experimentally deter- 
mined Co-0 bond lengths using the bond valence sum 
(BVS) method allows us to estimate the difference in 
the valance of the two Co-ions. For the calculation of 
the BVS we used b =0.37 and R =1.70 Hp!]. The 
BVS as a function of temperature for the two Co sites is 
shown in fig. [3^c). For high temperatures we find that 
the BVS shows some scatter that reflects changes in the 
mobility and ordering of Na atoms between C0O2 sheets 
(see below) and possibly to geometrical differences that 
arises from the Na-ion potential imposed on the C0O2 
sheet over the same temperature range. [IBJ However, be- 
low ~150K we find a small but measurable and consistent 
difference of e ~0.06e in the BVS for the two Co atoms. 
Although the difference is comparably smaller than what 
is found in conventional charge ordered systems, the sep- 
aration of the data into two values (one low and one high) 
below 150K is statistically significant. These data would 
suggest that below Tco=150K there is a separation of 
charge into Co 3 5+e and Co 3 ' 5-6 stripes running along the 
b— axis as shown in fig. [l](a). This charge ordered struc- 
ture is in agreement with the magnetic neutron diffrac- 
tion measurements, where the magnetically active Co- 
site would correspond to the Co(l) site with the slightly 
higher BVS and octahedral distortion. The values of the 
BVS obtained are in good agreement with both theoreti- 
cal predictions [6| and recently reported values at 10 and 
300K respectively (l5| which are also plotted on fig. [3^c) 
for comparison. The weak charge ordering found here is 
consistent with the relatively low resistivity of this mate- 
rial at low temperatures ( ~ 100 raVL cm at 2K)0, 0, @]. 
That the average BVS is approximately 3.3e reflects the 
mixed valent nature of this compound. [26| . 



IV. NA ORDERING AND BEHAVIOR OF NA-O 
BOND LENGTHS 

We now turn our attention to the behavior of the Na- 
layer for this composition. In fig. [4] we plot the tempera- 
ture dependence of the atomic displacement parameters 
Uiso (Debye- Waller factor) determined from our Rietveld 
analysis. We find that the Ui SO values of the O- and Co- 
atoms to be in general of the expected amplitude and 
show a linear behavior with temperature. The behavior 
of Ui SO for the Na-ions however is unusual in that there 
is a clear change in slope at 300K separating a low and 
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FIG. 4: (Color online) Isotropic atomic displacement 
parameters (Ui ao ) for the Na, Co and O-atoms determined 
from the Rietveld analysis of the NPD data. In the analysis 
the following constraints were used L/j so (Co(l))=i7i SO (Co(2)), 
( 7 iso (Na(l))={7 iso (Na(2)),?7 iso _(0(l))=[/ jso (0(2))=[7 jso (0(3)). 
Lines through the data are guides to the eye. A slope change 
in the t7j aa (Na) is evident around 300K. In the inset we 
show the temperature dependence of the (111) reflection 
in the orthorhombic Pnmm setting. This reflection is a 
superlattice reflections with respect to the parent P63 /mmc 
crystal structure and arises from the ordering of Na ions. 
The temperature T s ~460K marks the decomposition of the 
sample in hexagonal Na I Co02 and C02O3. 



a high temperature behavior, while for T>300K the Ui SO 
values for Na become large. Such behavior is indicative of 
a dynamical transition occurring at 300K involving only 
the motion of Na-ions, as similar signatures are absent 
for the Co- and 0-atoms[18]. Indeed such large vales of 
Uiso suggest that Na-ions may become mobile between 
C0O2 layers above 300K. 

For higher temperatures we find that our sample de- 
composes at T s ^460K to Naj;Co02 and C02O3. This 
transition is quantified by tracking the intensity of the 
(111) reflection as shown in the inset of fig. |U This re- 
flection is a superlattice reflection with respect to the 
parent P63/771771C structure and arises from the ordering 
of Na-ions [|[. Our neutron powder data measured at 
475K indicate the loss of this and other superstructure 
reflections and a return to VQ^/mmc symmetry with the 
addition of C02O3 reflections. 

Within this perspective we now look more closely to 
the temperature dependence of the Na-0 bond lengths 
show in fig. (Ha,b). For both Na sites the Na-0 bond 
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FIG. 5: Temperature dependence of and Na-0 bond lengths 
computed from Rietveld analysis of NPD data measured on 
the E9 diffractometer. Dashed lines are guides to the eye. 



lengths show a set of short bonds (~2.36A) and a set of 
long bonds (~2.44A). While the high temperature behav- 
ior is complicated by the high Na-ion motion as discussed 
above, on cooling below 300K the long bonds decrease, 
while the short bonds show an increase down to 150K. 
This correlated behavior in general indicates a displace- 
ment of Na-ions along the a— axis. At 150K we find that 
the small charge disproportionation in the C0O2 layer is 
not reflected in the Na-0 bonds. Surprisingly however we 
find that at T m \ a decrease of ~0.0lA of the long Na(l)- 
O bond and a correlated increase in the short Na(l)-0 
bond, while for the Na(2)-0 bond lengths an opposite 
and less clear effect can be seen in the data. 

We interpret these changes in the Na-0 bond lengths 
as reflecting changes in the electronic state of the C0O2 
layer at T m \. The nature of any coupling between 
changes in the electronic state of the C0O2 and Na can 
arises from the orbital configuration of the Co-ion it- 
self. It is argued by Kroll et al. [IH that the edge-sharing 
Co 4+ 06 octahedra are compressed along the c— axis re- 
duces the point group symmetry to D^d- The t 2g orbital 
of the Co 4+ is split in D 3 d as t 2g = a'l g + e' g , giving 
a fully occupied e' g and a half filed a'l g . The latter or- 
bital looks like a 3d z 2_ r 2 orbital and points along the 
c— axis|l9|. For the case of the higher valent Co(l)-ion 
(nominally Co 4+ ), this orbital would point in between 
the O-atoms and towards the Na(l)-atom as shown in 
the inset of fig-EJa). Since the charge disproportionation 
here is small each Co ions will have a similar electronic 
and orbital configuration. 

Therefore the orbital configuration of the Co-ions pro- 
vides a means to couple electronically the C0O2 and Na 
layers. The nature of the coupling is electrostatic and 
would arise from the occupation of the e' g orbital. We 
would expect that changes in the electronic configura- 
tion of the C0O2 sheet to be reflected also in the relative 



positions of the Na-ions as indicated by the Na-0 bond 
lengths. More precisely changes in the electronic state 
of Co should be clearest for the Na(l)-0 bonds as the 
Na(l)-ion sits directly above (or below) a Co(l)-ion. The 
same argument would suggest a less pronounced effect for 
the Na(2)-0 bond lengths as the Na(2) ion resides above 
and between C0O6 octahedra. This is indeed reflected in 
the bond length data where a strong response is found 
in the Na(l)-0 bonds and a less clear responce in the 
Na(2)-0 bonds [27|. At T m2 we find no clear evidence of 
changes in the Na-0 bond lengths. This is expected as 
the changes in charge separation at this lower transition 
are computed to be much smaller than those at T m iQ. 



V. ANOMALOUS BEHAVIOR OF LATTICE 
CONSTANTS 

In fig. [6ja-b) we show the temperature dependence 
of the lattice constants determined from Rietveld refine- 
ment of the NPD data. These data show a positive ther- 
mal expansion (TE) for the c— axis between 2 and 450K, 
but for the a— and b— axis we find an almost constant TE 
between 2 and 300K; here for T <300K linear TE expan- 
sion coefficients are — 9(3)x~ 7 /K and 1.1(4) x~ 6 /K for a 
and b respectively. Such small TE was also discussed in 
ref. EH 

for a much more limited number of temperatures 
and smaller range in temperature. For T > 300K there 
is a return to positive TE for both in-plane parameters. 
This crossover coincides with the dynamic transition in 
the motion of the Na-ions as indicated by the behavior 
of Uiso for Na. 

The behavior of the lattice constants for this x =0.46 
sample is in sharp contrast to our x =0.75 sample (shown 
in the inset fig. , were we find a positive TE for both a— 
and c— axes between 5 and 300K. Assuming that TE is 
dominated by acoustic phonons below 300K for Naj;Co02 
materials and whose frequency is relatively invariant be- 
tween x =0.75 to x =0.46, the T— dependence of the 
lattice constants for the x =0.75 sample (see inset in 
fig-E^f)) can be used to quantify the difference in the in- 
plane TE between these two samples. Here we define the 
term a = (a' — a^ e2; )/a / where a' = (y/3a + 6/2)/2, a and 
b are the orthorhombic lattice constants of the x =0.46 
compound and a/ lex is the hexagonal lattice constant of 
the x = 0.75 compound. 28] For the higher temperature 
data (T>300K) ahex was assumed to vary linearly with 
temperature. Here a represents an additional temper- 
ature dependent contribution to the expected in-plane 
lattice constants (as defined by ahex) and its tempera- 
ture dependence in shown is figEJc). For T> 300K a 
is near zero, although the data in this region are more 
limited. However for T<300K a increases with decreas- 
ing temperature and reaches a value of ~ 40 x 10 _3 % 
at 1.8K. The maximum value of a is ~ 50 x 10 _3 % at 
150K, close to T co . The decrease below this tempera- 
ture is due to an increase in the ahex for the x =0.75 
as seen in the inset of fig. HJc). We speculate that the 



6 




1 50 300 



100 200 300 400 
Temperature (K) 

FIG. 6: (Color online) (a,b) Lattice constants and unit cell 
volume (inset in panel (a)) obtained from the Rietveld refine- 
ment of NPD data measured on HRPD at ISIS (filled circles) 
and E9 (open circles) over the temperature range of 2 to 450K. 
In panel (a) and inset the solid lines represent the thermal 
expansion obtained from a fit to the data using the second 
order Griineisen approximation, (c) Temperature dependence 
of the additional component to the thermal expansion a for 
Nao.46Co02 obtained by subtracting the temperature behav- 
ior of the in-plane lattice constant of the x —0.75 compound. 
In the inset of panel (c) we show lattice parameters measured 
as a function of temperature from a hexagonal x =0.75 sam- 
ple. 



anomalous TE below 300K may be driven by electronic 
correlations induced by Na-ions. 



VI. DISCUSSION 

The measurement we present here suggest a picture of 
incipient charge ordering for near half doped Na x Co02. 
At high temperatures the ordering of Na-ions defines two 
different CoOg octahedra, one that is relatively undis- 
torted and one that more distorted. On cooling the 
differences in octahedral distortions between these two 
different Co-sites becomes larger and may reflect an in- 
creasing influence of the Na-ion potential on the C0O2 
sheet, as Na ions become more localized around their 



mean crystallographic positions. Indeed the near zero 
TE for the in-plane lattice constants coincides with a dy- 
namical transition in the displacement parameter Ui SO 
of the Na-ion. Further the physical meaning of a can 
be interpreted as a measure of the Na-induced electronic 
correlations onto the C0O2 layer which in this view sat- 
urate at T C o=150K. 

While our NPD work can correlate the distortion of the 
C0O6 octahedra and the Na-ordering even at high tem- 
peratures, it is not until Tco that we find evidence for a 
weak charge separation into stripes. Indeed it is possible 
that this pattern of charge ordering is present in the lat- 
tice from the onset as a direct result of the Na-ordering 
but it in effect is hidden by the dynamic behavior of 
the Na-ions. Therefore we argue that the charge order- 
ing emerges at low temperatures as Na-motion becomes 
more confined. 

More critically in terms of the physics of these ma- 
terials we demonstrate that charge ordering occurs at a 
higher temperature that the magnetic ordering and elec- 
tronic transitions at T m \ and T m 2 respectively. This is 
consistent with recent theoretical models that suggest 
that the Na potential imposes a degree of charge or- 
dering to the lattice [1, 0, 13] ■ Indeed Zhou and Wang 
[(| , suggest that as much as half of the expected charge 
disproportionation would occur at a temperature above 
TLi [291] . consistent with our observations. 

At lower temperature the changes in the C0O2 layer 
may be inferred indirectly by monitoring the Na(I)-0 
bonds. Here the orbital configuration of the Co pro- 
vides for charge density pointing directly to the Nat- 
ions thus providing a sensitive parameter to electronic 
changes in the C0O2 layer. Indeed changes in the Na- 
O bond lengths may be more sensitive than changes in 
Co-0 bonds as e g axial orbitals are empty. Our mea- 
surements find that that changes in the Na(l)-0 bond 
lengths correlate with the magnetic transition at T m \ 
suggestive of further changes in the electronic state of 
the C0O2 layer. It is predicted that charge separation 
is enhanced gradually below T ml [|| but the changes here 
are overall again small and may fall outside the limits of 
our sensitivity. At T m i we find no evidence of changes 
in the lattice or changes in the lattice symmetry. The 
prediction of a modulation of the amplitude of antifcrro- 
magnetically coupled spins as well as the charge within a 
Q +3.5+e gt r jp e i s much smaller than our detection limit 
(~0.02e)@. 

The structural observations at Tco we report here cor- 
relate with features in the charge dynamics. For example 
Quian et al. [201 ] report from ARPES measurements that 
with increasing T from the insulating region (were a clear 
gap is found) the size of the gap and the spectral weight 
around the gap decrease. Although the gap closes at T m \ 
the spectral weight does not completely vanish until ap- 
proximately I20K, a behavior that is attributed to the 
formation of quasiparticles that gain significant weight 
due to coupling along the c— axis. For similar temper- 
atures optical spectroscopy measurements find a broad 
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feature that is associated with fluctuating charge order- 
ing or a CDW in both anhydrous [2li [22| and hydrated su- 
perconducting samples [23]. These observation together 
with our structural measurements point towards a pic- 
ture where at 150K an incipient charge ordering forms. 

VII. SUMMARY 

In summary this work establishes that (a) the Na- 
ordering on the C0O2 layer is reflected in the octahedral 
distortion of the two crystallographically in-equivalent 
Co-sites and is evident even at high temperatures; (b) 
The charge ordering occurs below Tco=150K, a tem- 
perature higher than the magnetic ordering found at 
T m i=88K, consistent with theoretical models that sug- 



gest that the Na potential imposes a degree of charge 
ordering to the lattice [f| 0, E3]; ( c ) Below T co we find a 
weak charge ordering into stripes of Co 3 ' 5+e and Co 3 ' 5-6 
with a e ^0.06 e, a value in good agreement with that 
obtained from a Hubbard model using the Gutzwillcr 
approximation 0; (d) A dynamic transition in the mo- 
tion of Na-ions occurs at 300K and coincides with the 
onset of a near zero thermal expansion for the in-plane 
lattice constants of our Nao.46Co02 sample. 
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